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Abstract 

A crucial requirement for the standard interpretation of Monte Carlo simulations of sim- 
plicial quantum gravity is the existence of an exponential bound that makes the partition 
function well-defined. We present numerical data favoring the existence of an exponential 
bound, and we argue that the more limited data sets on which recently opposing claims 
were based are also consistent with the existence of an exponential bound. 



1 Introduction 

Theories that we hope could be connected to Euchdean quantum gravity can be inves- 
tigated through Monte Carlo simulations of discrete models, one of the most popular 
approaches being dynamical triangulations of simplicial quantum gravity. Here we con- 
sider the four dimensional theory without matter but with a cosmological constant for 
fixed topology 7V4 = 5^ |]l|, |2|, |3|, §]. The most attractive feature of this model is that there 
is evidence for a second order phase transition and therefore a continuum limit may exist. 
In this language the Einstein-Hilbert action is very simple. It depends only on the 
total number Ni of i-simplices contained in the triangulation T: 

S[T] = kiNi[T] - k2N2[T] . (1) 

The coupling constants ^^4 and k2 are directly related to the cosmological constant A and 
Newton's constant G, 

h = X-^^, k2 = ^, (2) 

G aG 

where a = arccos(j). The grand canonical partition function of the theory is 

Z{k^,k2) = 5]e-^[^] (3) 

T 

^ y^g-fc4Af4 y^ gfc2A^2[T] 

N4 T:\T\=N4 



where we have split the sum over all triangulations of A4 into a sum over all possible 
volumes (equal to the number of 4-simplices A'4) and a sum over all triangulations T 
with volume |r| equal to A^4, i.e. the second sum gives the canonical partition function 
Z{N4,k2). 

The question that we want to address in this note, and which has been prompted by 
Q and further discussed in |^, ^, is whether the partition function of such a model is 
actually well-defined, i.e. finite. Suppose that there exists an exponential bound for the 
canonical partition function, 

Z(iV4, A:2)~ 6^=4(^2)^4 ^ (4) 

for large N4 and some constant k2^{k2)- Then the partition function Z{k4, /C2) is finite for 
^4 > ^(^2) and divergent for k^ < ^4(^2). 

The question of the existence of an exponential bound for the canonical partition 
function is directly related to the asymptotic behavior of the number of triangulations for 
a given volume, J\f{N4), which might grow as fast as (5A'4)!. Since 2N4 < N2 < 4:N4, 

A/'(iV4) < Z{N4, /C2) < e^'=2A^W(A^4) if k2 > 0, (5) 

Z{N4, k2) < e^'"'^^JV{N4) < MiN^) if k2 < 0. (6) 

Hence the existence of an exponential bound on J\f{N4) implies the same for the canon- 
ical partition function for arbitrary k2, and if there exists an exponential bound on the 
canonical partition for a single value k2 >0 then it exists for all k2- 

In the next section we discuss how the absence or presence of an exponential bound 
manifests itself in the numerical simulations, but first let us summarize the history of the 
subject. Until recently the existence of an exponential bound for the number of trian- 
gulations for large volumes was considered very probable. The bound can be rigorously 



proven in two dimensions [^ . In three dimensions there is strong numerical evidence for its 
existence Q , and in four dimensions the numerical data still seemed reasonably consistent 
with that claim (e.g. Q, |, |, |). 

However recently the authors of H have claimed that their more detailed examination 
of new numerical data shows the absence of an exponential bound. (In ||5[ the coupling 
constants are kq = 2A;2 and K4 = A;4 — 2/c2-) Their main points are: 

1. In 2d everything is fine, and the bound is manifest. 

2. In 3d k^{N4) has strong power corrections to a (to be proved) asymptotic constant 
behavior. One does not see the asymptotic behavior, but the transient behavior is 
power like. 

3. In 4d A;|(A'4) is not constant, but changes far less than in 3d. The authors of |^] 
assume this means that in this case there is a logarithmic divergence of fel (while 
the behavior in 3d, far more abrupt, was assumed to be convergent) that implies a 
violation of the exponential bound. 

4. The situation in d ^ 4 is mainly discussed for k2 = 0, while for 4(i data is presented 
also for k2 = 0.25 and 0.50, all lower than /cg ~ 1.1, based on the idea that proving 
existence of the exponential bound at one value of ^2 ^ is sufficient for proving 
that the partition function is convergent everywhere. 

After [^, two papers have discussed the issue in further detail, both observing that a 
logarithmic and a small power fit are reasonably consistent with their data. The authors of 
ref. Q lean towards the logarithmic scaling, and they propose for this case an interesting, 
potentially still well-defined scenario, on which we will comment below. Ref. Q argues in 
favor of a power law approach to an exponential bound, proposing the ansatz of a leading 
power a = -g. They fit the old and some new data that look quite consistent with the 
ansatz. 

Here we present data for d = 4 and k2 = for the largest system size yet, and 
conclude that the extension to larger systems allows one to decide that there does exist an 
exponential bound — if one is willing to accept as a compelling evidence the fact that the 
ratio of x^ of the factorial fit to x^ of the exponential fit is 10. We discuss the strength 
and the implications of such a result. 

Points (1) through (4) remind us that there is an important issue to be clarified, but 
also make the numerically oriented physicist quite perplex. It is well known indeed that 
measuring logarithmic corrections is very difficult (for a classical study see for example 
pop). So in the case of point (||) it is impossible to exclude that there could be a logarithmic 
divergence, underlying a transient behavior dominated by power law corrections. The main 
point in all the following will be that we are discussing a transient region, with an overlap 
of different possible corrections, which could all be present at the same time. 

The same kind of perplexities arise about point (|3|). When numerical data taken in 
some limited range of the parameters are compatible with a logarithmic fit they can also 
be fitted by a power law with a small exponent. What one can achieve is putting a bound 
on the allowed value of the power, and in our case we will show that indeed a convergent 
power dependence has to be preferred over a logarithmic scaling. 

We will also discuss in more detail the different k2 regimes (the crumpled and the 
elongated phases), but already the evidence contained in Q makes the appreciation of 
point (^) quite delicate. Indeed, figure 2 of [^ (also Q, and later figure 2 of |^, to be 
compared) shows that in the two different phases k^(N4) is behaving in a very different 
way. In the elongated phase k^{N4) does not depend at all on the volume, while the 
residual dependence is all in the crumpled phase. It is possible that because of strong 



finite size effects in the crumpled pliase, reliable numerical data about the large volume 
limit can only be obtained above the apparent phase transition for k2 > k2 ~ l-l- 



Finally, in [11 1 hints are given toward the fact that in Ad the number of triangulations 



could be exponentially bounded. In three dimensions, there does not yet exist a proof 



(e.g. |12|). Independent of any analytical arguments, however (and, as far as we know and 
understand, we cannot yet rely on a rigorous proof), one should arrive at the best possible 
understanding of the numerical data. This point we will address below, stressing again 
that we are talking here about the nature of a transient region, where different corrections 
of unknown form may conspire to make the picture difficult to disentangle. 

2 Results for k2 = 

When simulating the system described by the partition function Z{ki,k2) (^), with variable 
volume, one finds that there exists a line k1{k2) in the plane of the coupling constants such 
that for any k2, if A;4 is larger than A;| the volume tends towards zero (towards six to be 
precise) , and if A;4 is smaller than kl the volume goes to infinity. The larger the deviation 
from kl the faster is the trend. For these reasons ^4(^2) is often called the critical line 
(although it has nothing to do with statistical criticality |Q]). This critical line at finite 
volumes represents the behavior one expects from the discussion of the exponential bound 
in the introduction. 

A typical Monte Carlo simulation is performed close to a fixed volume A'^4 and for fixed 
^2, while kii is kept close to k1 which is determined dynamically during the simulation. One 
then looks for a phase transition determined by some /c2 = ^2 i^ the limit of large volume. 
How k1{k2) is measured is a purely technical question, and the particular algorithm we 
use is described in detail in |^]. 

In figure 1 we show k1 versus ln(A'^4) for A;2 = (1/G = 0). Here we are mainly 
discussing the point /c2 = because this is where in |^] the strongest argument against 
an exponential bound is made. If instead of an exponential bound only a factorial bound 
holds, then one expects 

kl{Ni) = a + h\nN^ . (7) 

If an exponential bound e""^"^ to the canonical partition function exists, then 

kl{Ni) = a + hN^'' , (8) 

where N^'^ represents a natural polynomial correction to the exponential. One can argue 
for a = 1/4 [^, which allows a very nice fit to the data. Since there is not enough data 
to determine a reliably, setting a = 1/4 serves well enough to distinguish the exponential 
from the factorial fit. 

In [0|, a straight line (we draw our best logarithmic fit with dashes in figure 1) has 
been chosen as the best fit to data points between A'4 = Ik and 32k, corresponding to 
the absence of an exponential bound due to factorial growth. In [^, a polynomial fit for 
a = 1/4 is preferred for A'4 from 4k to 64k, representing the existence of an exponential 
bound which is approached only for still larger volumes. Volumes of Ik and 2k 4-simplices 
have not been included in the fits since they are likely to suffer from strong finite size 
effects. The data of [^, |^, suggest that we require more data at larger volumes for a 
reasonable estimate of the asymptotic behavior. 

It is just a matter of computer time to take data at larger volumes — more so than 
a matter of computer memory, and we estimate that our implementation uses only about 
one fifth the memory of the implementation of |0]. We have been able to get reliable 
data up to a volume of 128k (50000 sweeps, counting A'4 moves actually performed, for 
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Figure 1: A;| versus ln(A^4) for /c2 = 0. Values for A^4 are 4000, 8000, 16000, 32000, 64000 
and 128000. With the dashed hue we give our best logarithmic fit, with the solid line we 
give our best fit to a converging power, with a = .25. Both fits have two free parameters. 
The x^ of ths power fit is ten times better than the one of the logarithmic fit. 



A'4 = 4:k through 32k, about 30000 sweeps for A'4 = 64A; and 128k, 10000 sweeps are dis- 
carded for thermalization, the largest volume taking four months on a (shared) IBM/RISC 
workstation, autocorrelation time was on the order of 50 sweeps). 

It is remarkable that when superimposing our new points to the fits of ref. |M] they fall 
very well on the power fit (quite far indeed from the logarithmic divergence prediction) 
obtained from smaller triangulations. 

We have fitted our data for k2 = with the two forms (j^) and (P) , by setting the power 
a = i. They are both two parameter fits. Figure 1 is quite eloquent about the success of 
the two fits. The result is 

j^cilog) ^ 0.864 +.0277 In iV4 , (9) 

^c(po«,er) ^ ^252 - 1.3177V4"^ . (10) 

The power fit has a value of the x^ which is ten times better than the logarithmic one. We 
have also tried 3 parameter fits. In the power fit we have left the power as a free parameter, 
while in the logarithmic fit we have added a volume scale term N^, as in ln(A'4 — N^). 
Both fits improve quite a lot, but the power fit stays far superior to the logarithmic fit 
(the x^ ratio is now 3). While such a power fit (where the best power is now .36 it .04) 
matches perfectly the data points, the logarithmic fit is still not totally congruent to the 
data (we get N^ of order 3000, that is a reasonable scale for the transient behavior). We 
are not very confident in playing with many parameters, since the allowed corrections are 
of many different functional forms, and it is clear that with 6 data points they cannot be 
distinguished. We just take the results of the 3 parameter fits as further evidence that the 
power fit is superior to the logarithmic fit. Let us also note that indeed the best preferred 
power is surely not too small. 

The conclusion we draw is that the fits of the numerical data largely favor the existence 
of an exponential bound at k2 = over the presence of a factorial bound. 

3 Discussion 

What about the consistency of the numerical data? The first observation about the data 
in figure 1 should really be that there is a remarkable agreement in the data from four 
independent computer implementations considering that the underlying algorithms are 
somewhat similar but not identical. In fact, notice that even the data from [^ that lead 
to the claim about the absence of an exponential bound curves away from a straight line 
in the same way the other data sets do. 

Having analyzed in detail the situation for k2 = 0, we now turn to generic values of the 
coupling /c2- In theory, the existence of an exponential bound for any one value of A;2 > 
implies existence for all the others. But as is well known, but has not been discussed in 
detail in this context, there is an important practical difference between the phases for k2 
below and above the critical value ^2 ~ 1.1. For large positive k2 the simplicial complex 
is in an elongated phase with an intrinsic dimension close to two, while for negative k2 
the intrinsic dimension diverges to infinity and the simplicial complex becomes extremely 
crumpled. One of the most intriguing and attractive features of simplicial quantum gravity 
is that at /c2 the intrinsic dimension is close to four |l|, 13 1 (for simplicity we ignore here 



the problem of giving the best definition of the intrinsic dimensionality of the system). 

The point is that the two phases are not only different, but there is a genuine asym- 
metry. Note that at /c2 the intrinsic system size for A'4 = 10, 000 is of the order of 
(10,000)^/^ = 10, while at A;2 = it is (10,000)^/^° ^ 2.5. Therefore, what constitutes 
a large volume that guarantees the absence of finite size effects depends very sensitively 



on the value of k2 |^]- For example, the asymmetry in the susceptibility present in these 
systems may be due to such effects. 

With regard to the discussion of the exponential bound one should therefore consider 
the whole k2 range. Such data already exists in [Q, Q and were improved upon near the 
transition in |^] but were not considered in p, 0]. For concreteness we show in figure 2 a 
plot of X'^{k2) versus Aq ~ l/G for A'^4 = 4/c, 8k, 16k based on ||3|, which for our purpose is 
better suited than the more accurate data of @ since figure 2 extends to extreme values 
of k2- The constants are defined by the relations 

k2 = 2ttXo , k4 = X + lOaXo . (11) 

There is a definite volume dependence for k2 < k2 while above the transition no volume 
effect is discernible. The linear transformation from k2 and /C4 to the cosmological constant 
A is useful for magnifying the volume dependence which is invisible in this range of coupling 
constants for k2^[k2) [^. This is discussed [^, but even when explicitly looking for a small 
volume dependence for k2 clearly above /c2, none is found. In this region the plot analogous 
to figure 1 appears to be a perfectly horizontal straight line, i.e. there are no detectable 
polynomial corrections to the exponential bound. 

The discussion can be taken one step further by noticing that the critical value k2 of 
k2 moves to larger values with increasing volumes |j|, ^. For larger volumes at /c2 = 
the finite size effects become even more pronounced (internal dimension up to 50). Given 
that for extreme values of k2 the simplicial complex freezes and ^4(^2) becomes a perfect 
straight line with different slopes, the shift in ^2 keeping the part k2 > fcj in figure 2 fixed 
translates directly into (part of) the volume dependence in the range k2 < k2. When 
this effect constitutes the significant part of the volume dependence (for large enough 
volume), then the volume dependence of the critical value /c2 can be estimated by the 
volume dependence of kl for a small enough but fixed value of k2- In particular, if there 
is no exponential bound, then ^2 ^ 00 with A'4 — > 00. 

It is instructive to examine the condition for the critical line in the Monte Carlo 
simulations (here we follow |^, Q). Consider the ergodic random walk in the space of 
triangulations of S'^ consisting of the five standard moves, where for the move of type i 
an z-simplex is replaced by a (4 — z)-simplex. On the critical line, the average volume A'^4 
is constant, and therefore N2 must also be constant since it is bounded. This means that 
the average variations 6Nj must vanish, 

4 

5A^~5]AAr,(i)p, = 0, (12) 

1=0 

where pi is the probability with which a move of type i is performed on average, and 
ANj{i) is the change in Nj due to that move. Since the moves are independent, we obtain 

P0=P4, Pl=P3, (13) 

on the critical line. 

Since the action is linear in N2 and A'4, and since the moves are local, we can be more 
specific about the conditions on the pi. The pi can be chosen to be 

K = [e-^^«]pr. (14) 

The bracket is the Metropolis weight, its key feature being that it depends only on the type 
of move and not on the Ni or the triangulation in general. While the action looks quite 
trivial, all the non-trivialities are hidden in the probability pf^" for a move to be allowed 
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Figure 2: A'^ versus Aq. Indications of a phase transition are found near Aq = 0.18. 



by the geometric constraints on the triangulation. (Detailed balance is incorporated in 
the way the moves are chosen. Potentially, there is a factor of order 0{1/N4).) 



Therefore (13) is equivalent to 



kl = h2-lnpr, (15) 

geo 

kl = 2k2-\n%ro, (16) 

where we have used that p^^° = 1. The question of the existence of an exponential bound 
has therefore been translated into the question whether there exist appropriate bounds on 
the Pi = Pi{ki, A:2, A'4) which are independent of A'4. 

First of all, po < 1 implies that k^ is bounded from below by 2.5k2- The hard part is to 
find a suitable lower bound on po, for example, and although it may be possible to do so 
by some more detailed analysis of the space of triangulations, we do not have a conclusive 
argument. Notice that since moves of type 4 are always allowed, we have that pq > 0. 
However, a naive counting of possible moves of type and 4 around a fixed background 
triangulation gives po ~ l/A'4, which would be the divergent scenario, but the same kind 
of counting would also make 2d divergent. The counting is, of course, difficult because 
moves of type 1, 2, and 3 may change the geometric constraints. 

Coming from the numerical side, it is quite suggestive that e.g. the data for A'4 = 4k 
in figure 2 corresponds to kl{k2 > 4.0) = 2.497/c2 + c and kl{k2 < -4.0) = 2.002A;2 + d. 
This means that in the extreme ^2 regions the relevant geometric probabilities must be 
independent of k2 (combining (15) with ( p^ ) gives a factor of exp(A;2/2) for the opposite 



side). 

Considering the general structure of the phase diagram, the volume dependence can 
also be understood on the level of the random walk as follows. It is the moves of type 
4 that drive the system into the elongated phase (AA''2(4)/AA'4(4) = 2.5), while moves 
of type 1 drive to the crumpled phase (AA''2(3)/AA'4(3) = 2.0). Depending on /c2 the 
random walk is driven towards one of the bounds in 2 < N2/N4 < 4. One of the two 
possible phases, the elongated phase, is therefore characterized by low order vertices, and 
the average order does not depend on A'^4 since a maximal elongation can be obtained for a 
rather small number of simplices. Hence Pq^°, which is the ratio of the number of vertices 
of minimal order to the number of all vertices, is expected to be independent of A'4 in the 
elongated phase. On the other hand, in the crumpled region the average order of vertices 
is driven towards large values, and the average order will grow with A'^4. Hence Pq'", which 
is defined by the low order tail of the vertex order distribution, goes to zero with N4 in 
the crumpled phase. Equation (|l5|) gives the corresponding volume dependence of k^. 

In conclusion, when looking for evidence for an exponential bound in the numerical 
data of simplicial quantum gravity in four dimensions, one should take the whole range of 
k2 into account. If one insists on looking in the crumpled phase at ^2 = 0, the numerical 
data strongly support the validity of an exponential bound. 
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